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design  results  in  making  the  cost  of  ths  bridge  prohibitive  or  the  risk  of 
causing  flood  damage  excessive,: 

In  the  past „  studies  by  the  U0S,  Geological  Survey   and  the  Bureau 
of  Public  Roads  pertaining  to  backwater  effects  caused  by  constrictions  have 
considered  shapes  of  openings  sue  s  as  those  produced  by  s  r i igesc 
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However,  very  little  has  bean.  don<s  in  the  way  of  making  a  systematic  study 
of  the  hydraulics  of  river  flow  under  the  various  shapes  of  arch  brides « 
The  arch  is  unique  in  that  the  surface  width  within  the  barrel  of  the  arch 
decreases  with  a  corresponding  increase  in  stage0  The  purpose  of  t> 
search  is  therefore  to  study  the  hydraulics  of  arch  bridges  30  as  to  com- 
pensate for  the  loss  of  efficiency  at  high  flows,  and  to  provide  a  method 
for  computing  the  backwater  upstream  of  the  bridge0  In  addition,  a  prac- 
tical method  of'  making  indirect  measurements  of  flcod  flows  at  arch  bridges 
is  proposed*. 

A  project  was  initiated  in  the  Hydraulics  Laboratory  at  Purdue 
University  to  study  this  problem,,  It  is  sponsored  by  the  Indiana  St 
Highway  Department  in  cooperation  with  the  U«S0  Bureau  of  Public  Roads-. 

HISTORY: 


The  earliest  systematic  laboratory  investigation  of  flow  through 
contractions  in  open  channels  was  performed  by  £<.  Wo  Lane0A  He  related  the 
discharge  and  the  water  surface  elevation  through  to  contraction  by  means 
of  empirical  discharge  coefficients,  and  indicated  that  there  may  exist 
seme  relationship  between  these  cosffisients  and  the  ratio  of  the  max? 
backwater  depth  produced  by  the  contraction  to  the  normal  depth  of  r 
without  the  contractions  This  ratio  is  referred  to  as  tha  backwater  ratio-, 

In  1955.?  Kindsvater  and  Carter  presented  a  practical  solu' 
of  the  discharge  equation  by  an  extensive  experimental  investigation.. 
applying  correction  terms  for  various  geometric  conditions  "to  a  stan< 
discharge  soefficient,  the  method  can  be  applied  to  a  wide  variety  of  bound- 
ary conditions*  A  detailed  description  of  the  internal  and  external  flow 
characteristics  was  givene 


the  same  year.  Ho  J«  Tracy  and  R„  Wo  Carter^presented  a  com- 
panion paper  to  the  on  3  by  Kindsveter  and  Carter.,  In  it  they  gave  a  method 
of  computing  the  nominal  backwater  due  to  open-channel  tfOjintriofctoiVJ  a  The 
practical  solution  was  based  upon  empirical  discharge  coefficients  and  a 
laboratory  investigation  of  the  influence  of  channel  roughness,  chan 
shape,  and  constriction  geometry.  Their  study  ws  limited  to  single.  spans 
deck-type cjy£i8jti*itttian3  and  to  steady,  tranquil  flow,.  C.  ?»  Izzard,  *in  his 
discussion  of  this  paper,  pointed  out  that  the  backwater  ratio  is  definite 
a  function  of  the  normal  depth  Froude  number  at  the  constricted  section., 
Also  he  questioned  the  use  of  the  backwater  ratio  concept  when  the  head  loss 
between  the  section  of  maximum  backwater  and  the  vena  contracta  is  large 
compared  to  the  approach  velocity  heado 

The  combined  work  of  Kindsvater,  Carter P  and  Tracy  was  organized 

e 
into  a  UoSo  Geological  Survey  Circular/  which  presented  a  method  for  deter= 

mining  peak  discharges  at  abrupt  contractions 0  The  discharge  estimate  was 

to  be  made  from  a  survey  of  high-water  marks  and  channel  characteristics,. 

Although  the  method  applies  well  to  deck  type  bridges*  there  is  no  cirect 

application  for  using  the  method  when  an  arch  bridge  is  U3ed  to  make  an 

indirect  me  asurement o 

In  October,,  1957,  the  Colorado  State  University  in  cooperation 

with  the  UoSo  Bureau  of  Public  Roads  published  a  bulletin  by  Ho  KD  Li'V 

Jo  No  Bradley  and  £a   Jo  Plate  entitled  "Backwater  effects  of  Piers  and 

Abutments" o  A  rigorous  and  extensive  investigation  of  the  backwater  effects 

of  piers  and  abutments  has  been  givene  The  paper  includes  a  complete  analysis 

of  the  energy  losses  through  the  constriction  In  the  end3  an  approximate 

simple  method  of  analysis  is  provided  for  the  highway  engineer  to  use,,  The 

general  principle  of  the  method  is  the  conservation  of  energy 0  A  number 


of  graphs  based  upon  laboratory  data  were  developed  for  determining  the 
maximum  backv/ater  and  the  differential  level  of  water  surface  across  bhe 
embankment*  Thi3  method  was  reproduced  in  a  bulletin  published  by 
Bureau  of  public  Roads 'in  October  <,  195&,  Much  of  the  work  done  at  Colorado 
has  been  used  as  a  comparison  to  the  present  research  and  reference  to  it 
vrf.ll  be  made  throughout  this  texto 

Ho  Ro  Valient ine°re ports  on  tests  performed  to  study  the  char= 
acteristics  of  flow  in  a  rectangular  channel  with  symet ideally  placed,  sharp 
edged  constriction  plates  placed  normal  to  the  flow0  The  flow  is  related 
to  the  upstream  depth  by  means  of  a  weir  type  discharge  equation* 
experimental  coefficients  were  found  to  depend  upon  the  geometry  of  the 
constriction  arid  the  Froude  Number  of  the  unconstricted  flow*  The  condi- 
tions which  prodiice  an  increase  in  upstream  depth  were  investigated  and.  the 
extent  of  the  increase  evaluated 0 

Some  recent  work  done  at  Lehigh  University  te.IIs"about  the  effects 
of  placing  spur  dikes  on  the  upstream  side  of  a  bridge  contraction^,  These 
dikes  are  designed  to  increase  the  hydraulic  efficiency  of  the  bridge  cross- 
ing* The  paper  presents  a  good  qualitative  description  of  the  energy  loss 
through  the  contraction,. 

H\isain10carried  out  a  preliminary  investigation  upon  which  the 
present  research  is  based*  He  studied  both  two  and  three  dimensional  semi= 
circular  arch  openings  in   a  smooth  flume c  General  centerline  surface  pro- 
files were  obtained  and  recommendations  for  future  studies  were  made*  A 
dimensional  analysis  of  the  problem  was  presented,, 

Sootqr-'-developed,  for  the  two  dimensional  case,  both  exact  and 
approximate  solutions  of  the  discharge  equation,,  He  also  continued  the 


5o 

small  flume  tests  started  by  Husaino  Included  were  two  and  three  dimensional 
semicircular  and  two~dimensJ.onal  segment  tests  in  a  smooth  and  rough  channel0 
Several  curves  relating  the  backwater  ratio  to  the  normal  depth  Froude  Number 
were  presented  for  several  arch  diameters,,  A  more  detailed  presentation  of 
Mi«  i-eiuiia  •••£  these  tests  will  be  discussed  in  a  later  section,.  Much  of 
Mr*  Sooky's  data  has  been  reanalyzed  to  fit  more  recent  techniques,, 

ANALYSIS ; 

Figure  1  shows  a  definition  sketch  of  the  effect  of  a  char. 
constriction  on  the  water  surface  profile^  Section  views  B  and  G  illustrate 
the  two  types  of  centsrline  surface  profiles  obtained  with  mild  and  steep 
slopes »  The  most  generally  occuring  situation  which  appears  in  actual  prac= 
tice  is  idealized  by  section  B0  The  depth  y0  is  at  a  point  far  enough  up- 
stream such  that  the  flow  is  basical3.y  unaffected  by  the  M]_  backwater  curve-, 
yi   is  the  point  of  maximum  backwater c  y2  is  at  the  section  of  ndnimum  .jet 
area  or  the  vena  contracts,,  yo  is  the  minimum  water  depth  of  the  regain 
curve.,  and  y^  is  again  at  a  point  sufficiently  downstream  from  the  contrac= 
tion  v/here  the  flow  returns  to  the  normal  depth. 

For  any  physical  problem  such  as  this?  a  dimensional  analysis 

is  convenient  for  the  purpose  of  guidance  and  interpretation  of  a  test; 

program*  In  this  manner  the  basic  variaoj.es  can  bs  grouped  into  dimea- 
sionless  quantities  and  their  relationships  investigated.,  In  the  proa. 

at  hand,  it  is  desired  to  determine  the  maximum  water  depth  upstream  of 
the  constriction o  It  is  assumed  that  the  variables  which  govern  the 
water  superelevation  may  be  grouped  into  three  catagories  as  follows: 
(The  reader  i3  asked  to  refer  to  Figure  1  for  an  illustration  of  the  ter- 
minology,, ) 


a„)  For  the  fluid 

V»  the  absolute  viscosity 
£P  density  of  the  fluid 
gP  acceleration  of  gravity 

bo )  For  the  strearr.  flow 

y1?   maximum  water  depth  upstream  of  constriction, 

(section  1) 

yrt,  the  normal  depth  of  flow  in.  the  approach  channel, 

(section  o) 

Vo,  the  velocity  of  flow  at  normal  depth0 

n,  Manning's  roughness  coefficient  of  the  approach 
channel o 

/5ih,  the  maximum  water  surface  drop  across  the  cor:3tric= 
tion0 

Co)  For  the  constriction 

Aq — The  total  normal  depth  flow  areac 

Aa""=The  flow  representing  that  portion  of  the  A  that 
'  passes  through  the  bridge  without  contraction.. 

Hence,,  from  the  abo-i'e  list  of  variables 5 

y.  -  i  (y*>  y*  >  «>  ^j  ~*  >  e  >#  >  ^  k)    :) 

Buckingham's  theorunr states  thut  in  a  physical  problem  including 
n  quantities  in  which  there  are  m  dimensions,  the  ouantities  may  bs  arranged 

into  (n-m)  dimansionless  parameters 0  With  the  mass,  length  and  time  systems 
of  unit?-  the  n--m  or  seven  dimensionless  rr  parameters  are  as  follows  t 

Inverting  the  first  two  parameters 

In  equation  (3)  the  term  M>«$»is  equivalent  to  the  square  of  the 
normal  depth  Froude  Number-  Also \tyjf  is  the  Reynolds  Number o  It  is  well 
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Figure     I 


known  that  gravity  forces  are  predominant  in  open  channel  flow  wher 
viseous  forces  play  a  secondary  rule.  The  Reynolds  Wumber  may  the 
be  disregarded  for  determining  yj^o  Furthermore,  by  assuming  that  the 
shape  of  the  water  surface  downstream  does  not  affect  materially  the  shape 
of  the  water  surface  upstream,  the  term '~r  can  also  be  eliminated,,  By 


combining  the  ratios  4f  and  4*  into  4&  and  excluding  the  above  mentioned 

%  %*  A.$ 

terms  equation  (3)  becomes 

The  backwater  ratio  is  therefore  expected  to  be  a  functior 
the  normal  depth  Froude  Number,  the  channel  roughness  and  the  ratio  ^%/a  c 

The  channel-contraction  ratio  (m{ )  is  defined  in  the  present 
research  as  that  portion  of  the  total  normal  depth  flow  which  can  p? 
through  the  bridge  waterway  without  contractiono  By  definition  it  5.  s 
equivalent  to  the  ratio  /W//L  obtained  from  the  dimensional  analysis „ 
Along  with  the  normal  depth  Froude  Number s   the  contraction  ratio  is 
haps  the  most  critical  variable  in  the  problem,,  As  defined t   the  contrac= 
tion  ratio  is  a  hypothetical  term  which  owes  its  significance  only  to  the 
geonetry  of  the  constriction* 

Referring  to  figurs  2,  for  the  rectangular  case,  the  total 
flow  is  that  flow  in  area  ADEH,  and  the  flow  that  passes  through  the 
bridge  opening  without  contraction  is  that  represented  by  the  area  E 
Therefore  the  contraction  ratio  m6  is 

m'  =  <j/Q  (5) 

If  we  assume  that  there  is  a  constant  uniform  velocity  V0  across 
the  whole  normal  depth  section,  equation  (5)  becomes 

nV    -    VQ  •   A>V%V:  Va-    bV«y/   %       (6' 


JJawavor,  fcr  an  arch  bridge,,  as  shown  also  in  figo  2,  the  sur- 
face width  will  be  different  for  each  and  every  normal  depth  y0o  There- 
fore in  the  same  manner  ? 

m'  -  V«  "  AiV'/AQV,   "  Va,  (?> 

The  ratio  of  the  two  areas  is  clearly  not  equivalent  to  b/Bo  (For  sim- 
plicity, b/B  is  hereafter  defined  by  the  symbol  m„) 

For  that  portion  of  a.  serai-sircular  arch  with  radius  r  and  depth 


13 

yo  the  area  becomes 

Aj  -  f°i\/rz-  ya  <fy     -    Z  ify^-yo2    +    r'sirVty} 


(8) 


The  arch  shown  in  figure  2b  has  a  radius  r  and  springline  width 
b0  The  arch  has  been  superimposed  upon  flow  area  of  depth  y©o  The  center 
of  curvature  is  at  a  distance  d  below  the  springline  of  the  archo  The 
flow  area  (AQ)  of  the  rectangular  channel  is  By©*  while  the  area  repre- 
senting that  flow  through  the  arch  is  given  by 

\=  frf^fy  -  t2^1^6^  (9) 

Now9  equation  7  becomes 

m'  «  _/H_  =        DV/r^-D^^-  rasin"D/r   _    4\/rx-<lx  -f-  rW^/r  flQ) 
A$  By»  8y0 

By  algebraic  manipulation ,  eq«  (10)  can  be  reduced  to  a  form 
containing  several  dimensionless  ratiosc  The  result  of  this  reduction 
is 

m'»  m  Cm  (H) 


10c 


v.' n  ere 


&:ic. 


with 


and 


m=  b/B 


C  =  — 


r>  d/r 


Li  the  form  of  equation  (II)  the  value  of  ia»b/B  is  ad jus : 
for  the  particular  arch  by  an  amount  equivalent  to  Cm  such  that  m- 
the  same  as  the  ratio  of  AjtoA^o  Kindsvater,  Carter  &  Tracy* and  Liu? 

1  assfc>i5»s*3  lave  defined  the  contraction  rations  Amply  b/B  or  l=b/Bc 
In  the  more  general  case-,  equation  (11)  can  be  used  for  vertical  abut- 
ment bridge  piers  as  idealized  in  figure  2a.  \>y  using  a  value  of  Cm 
unityo  Also,  previous  writers  have  stated  that  the  contraction  ratio 
is  equivalent  to  a  ratio  of  the  conveyances  in  the  contracted  and  un- 
contracted  regionso  The  authors  feel  that.,  as  defined^  m9  is  more  truly 
a  ratio  of  areas  rather  than  conveyances,  since  the  conveyance  empi 
both  the  hydraulic  radius  and  a  roughness  coefficient „ 

In.  the  general  case,  the  values  of  of  and  Scan  take  on  numbers 
within  certain  limits,  before  the  normal  depth  will  submerge  the  ere 
of  the  arch.  The  limits  are  as  follows 
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DEFINITION     SKETCH     FOR     THE    DEVELOPMENT    OF 
THE    CONTRACTION     RATIO 
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Figure     3 


'when  6  =  0  9   the  case  of  a  semi=>circular  arch  with  the  center 
of  curvature  at  the  springline  existSo  When  J5=/  the  arch  doe3  not  exist  „ 

The  values  of  Cm  have  been  calculated  for  several  <x's   and 
and  are  summarized  in  the  graph  of  figure  3o  The  submergence  limit  re- 
presents the  upper  limits  of  both  <=*  and|3  o  The  segment  arch  which 
a  constant  radius  arch  with  its  depth  below  the  springline  (i»e„  ^>c  ) 
can  b<?.  used,  as  an  arch  in  its  own  right  or  as  an  approximation  to 
elliptical  or  a  multiple  radius  archo  The  value  of  m'  for  the  latter- 
two  cases  could  also  be  determined  directly  from  eq0  7°  However,  they 
have  not  been  worked  out  in  the  present  researcho 

An  approximate  solution  of  the  discharge  equation  in  a  rectan- 
gular channel  with  a  sharp  crested  semi-circular  constriction  was  obtained 
and  is  expressed  in  terms  of  an  infinite  series  of  powers  of  the  ratio  )}>La 
With  reference  to  figure  \   the  Bernoulli  theorum  gives: 


q~   fvdA   -      S^W^gfy.-h)   '  2|/rl-h*~dh 


(■U) 


Expanding  equation  (14;  into  a  series  and  integrating  term  by 
term  and  making  use  of  the  fact  that  2r«=bs 


«  *  c*fidT*  !J-h  b  f1  - °-im(ff '  0-O'77(f)'*  - 


This  may  be  written  as 
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Q  =   cy,wb  M 


where 


and 


il, 


c  •    Cd  f4]/'Z3 


M  = 


"  o.\2<M$f-  o.oqjf^f 


+  --■ 


(15) 


(16) 


(17) 


(18) 


The  discharge  in  a  rectangular  flume  may  also  be  expressed  by 

Q   -  VoAo   -  Fo^T  e  u3/2-  (19) 

•fhere  F0   =     V°  A/gfa 

is  the  Froude  Number  o.f  the  undisturbed  normal  depth  flowc     Equating  (.15) 
•and  (19)  and  solving  for  the  coefficient  of  discharge 


-  .\2ML- .     F.      /  V.  )"  (20) 


Since  m=  ir>!/'£t~eq-  20  ir.ay  be  rearranged  such  that  the  backwater  rati: 


becomes 


■J-z 


1 7  Cj  rn 


■--, 


(21) 


Typical  values  of  the  d  t  Ci  az*e  ohoi4n  ia 

figure  1?  which      the  results  the  two-dimensional  semicircular  arch 
tests  in  the  rough  rectangular  channel*  It  is  interesting  to  note  the  limit- 
ing  conditions  of  the  discharge  coefficient  as  a  ^goes  from  zero  to  ; 
For  a  two-dimensional  ideal  orifice,.  Street erioshows  that  the  application 
of  complex  variables  to  the  "Schwarz-Christoffel  Theorum"  (better  known 
as  the  theory  of  free  streamlines)  leads  to  an  ideal  discharge  coefficient 
of 


2.k> 
+  2b   "  "¥+T 


The  coefficient  of  discharge  curves  of  figure  15  seem  to  converge  to 
o6ll  shovdng  that  this  is  a  limiting  value  of  Cd  as  ms  approaches  zeroc 


15  o 
When  m*  is  a:ual  to  unity*  Cm-1  and  b/B-  lo  Therefore 
there  is  no  contraction  at  alio  If  there  is  no  contraction,  then  lJ 
and    M     «  /   „  Also  J2Jp-  =  -^.Si  which  is  approx»  unity*  There- 
fore eq0  21  becomes 
at  m9"! 

/V  ,    \*/3 

^  -      ("Vci)  <*> 

or        Cj  »  /F© 

This  can  also  be  seen  in  the  graphs  of  the  discharge  coeffi- 
cient vs  the  contraction  ratio  for  the  two-dimensional  case* 

It  has  been  observed  by  the  authors  that  the  equations  derived 
by  several  different  investigators  for  the  backwater  ratio  produce 
various  constriction  geometry's  seem  to  have  a  basic  similarity  c,  as  an 
example,,  equation(21)  in  the  present  test  for  yjij   appears  to  be  some  func«= 
tion  of  (fe/^y)o  An  equation  for  the  backwater  ratio  given  by  Valent:' 
for  lateral  con  triction  plates  is 


(24) 


/^  V  C  (i  -m)  /  WKm1  J 

wliere  p.  ~  k  u  i  /•-.  \ 

m  -   V1^  *      !  ~  q     *    l  ~  m  C^m  *  0 

Do  ' 

Also  Liu^presents  an  empirical  formula  for  a  two-dimensional  vertical- 
board  model 

(k  )' =  4-483  K*|> '  !-(2-5  -M)J  +  i       (25) 


where  M"b/B=ml>     (Cm»l) 
Considering  only  the  leading  terra  1/M  in  the  quantity  in  bracket^  ( 
b  scomes 


K  *  *f$")  (25a) 


It  appears  that  v;ith  the  proper  interpretation  of  the  variables, 
namely  m"  &  VF0  8   the  results  of  tests  performed  on  different  geometric 
shapes  of  bridge  openings  should  produce  the  same  results0  For  instance, 
a  vertical  abutment  deck~type  bridge  may  physically  appear  completely 
different  than  a  semi--circular  arch  bridge 0  However.,  hydraulicaily  3peak= 
ing  if  they  have  the  same  contraction  ratio  m" ,  they  should  produce  the 
same  backwater  ratio0  The  limitations  of  the  assumption  must  necessarily 
lay  in  the  fact  that  both  bridges  must  have  the  same  eccentricity,  skew- 
ness  and  entrance  conditions 0  It  is  believed  that  this  concept  applies 
equally  as  well  to  multiple  span  bridges „  An  sttampt  has  been  made 
compare  the  two-dimensional  semi-circular  test  results  of  the  author,  the 
segment  data  obtained  by  Sooky,  and  the  V3  data  as  given  by  Liue6  The 
results  of  this  comparison  will  be  shown  and  discussed  in  a  later  section^, 

I^ERI^IfTAL^T^UF: 

Ac  Main  Testing  Facilities 

For  the  purpose  of  preliminary  testing^  a  small  variable  slope 
flume  6"  wide  and  12"  long  was  built o  The  channel  sides  and  bottom  were 
constructed  of  lucite  and  carefully  alligned  by  means  of  adjusting  screws c 
The  slope  of  the  flume  was  contr  lied  by  a  hand  operated  scissor  jack  at 
the  lower  end  of  the  flume 0  An  aluminum  I-beam  mounted  horizontally  above 


the  flume  served  as  a  track  for  the  mechanical  and  electric  point  g> 
used  in  obtaining  the  water  surface  measurements 0  The  elect vie  point 
gage  consisted  of  two  metal  points  that  wei-e  hooked  up  to  a  set  of  I 
teries  ard  a  galvanometer 0  A'hen  the  second  metal  point  would  make  contact 
with  the  water  surface,  the  circuit  would  close  and  the  galvanometer 
deflect e  The  flow  was  metered  by  a  1  inch  orifice  plat      2  inch  supply 
line*  Two  and  three  dimensional  testa  were  run  in  both     oth  anc"  rough 
flume 0  For  the  rough  tests g   the  wails  were  lined  with  copper  wire  r 
of  16  meshes  per  inc 

The  majority  of  the  tests  reported  here  were  performed  in  * 
larger  2  foot  by  fj  foot  by  64  foot  all  steel  titling  flume 0  The  slope 
was  controlled  by  six  screw  jacks  that  v/ere  designed  and  installed  such 
that  the  rate  of  rise  and  fall  of  each  jack  per  turn  of  a  single  drive 
shaft  was  proportional  bo  the  distance  from  the  pivot  point  of  the  flume0 
The  jacks  were  driven  by  a  common  motor  and  gear  reducer,,  The  motor  was 
operated  by  a  raises  lower  and  stop  switcho  A  revolution  counter  way 
attached  at  one  end  of  the  drive  shaft  and  the  actual  slope  of  the  £ 
bed  was  related  to  the  number  of  revolution  and  tenths  of  revolutions 
of  the  sharto  In  this  manner  a  change  of  slope  with  an  accuracy  of 
£  O.OOOOOZS   feet/feet  was  easily  accomplished  in  a  matter  of  rainutes0 
At  the  discharge  end  of  the  flume  an  adjustable  sharp  crested  rectangular 
weir  made  of  lucite  was  installed0  A  catchment  box  was  made  to  eliminate 
any  splasho  The  box  discharged  directly  to  the  suiap0  An  8  foot  oy 
f  oot  head  box  was  equipped  with  an  elliptical  transition  to  provide  a 
smooth  change  as  the  water  flowed  into  the  flume.  The  ,-.ead  box  also 
contained  several  screens  and  one  large  stone  baffle t,  A  skiaming  board 


13o 
which  floated  on  the  water  was  installed  at  the  flume  entrance  in  order 
to  eliminate  surface  waves*  For  a  raore  complete  description  of  the  test- 

~i'«ii   flu1116  3(2®  reference  15  o 

The  x^ater  was  taken  from  a  large  recirculatory  sump*  One 
2000  GPM  puinp  and  one  300  GPW  pump  fed  the  head  hox0  The  actual  inf 
was  metered  by  2  venturi'So  A  complete  layout  of  the  £~±me   and  the  \.;ater 
supply  system  i3  t_bown  in  figure  4o 

Aa  aluminum  instrument  rack  was  mounted  on  adjustable  stain= 
less  sceel  guide  rails  running  the  length  of  the  flume,.  On  the  rack  was 
mounted  an  electric  point  gage  and  a  1/2  inch  Prandti  Tube,,  The  st 
of  the  point  gage  was  marked  in  millimeters  and  was  equipped  with  a  vernier 
which  read  to  a  tenth  of  a  millimeter,,  The  prandti  Tube  v;as  the  typo 
used  normally  for  air0  It  was  connected  to  an  inverted  U  manometer  which 
had  a  fluid  of  specific  gravity  o#10o  Two  surveyors  tapes  were  used  bo 
determine  the  location  of  a  particular  reading «  One  was  installed  length- 
wise  on  the  flume  wall,  and  another  transverselly  on  the  instrument  racko 
The  raut  along  with  the  point  gage  and  Prandti  Tube  3»n  be  seen  in  fi- 
gure £a0  In  addition  a  JjO-tube  manometer  stand  was  installed  to  obtain 
rapid  measurements  of  the  surface  geometry 0  Fifty  piezometer  taps  located 
at  points  along  the  centerline  and  1  ft„  and  2  ftc  right  and  left  of  the 
centerline  were  hooked  up  to  the  manometer  banke  The  bank  was  constructed 
so  that  it  could  be  tilted  to  a  45°angle  and  was  illuminated  from  the  insidec 
Figure  5b  shows  a  picture  of  the  c^pleted  manometer  stancL 

There  were  sixteen  models  used  in  the  testing  program,,  They 
were  designed  for  specific  values  of  b/B  and  L/b  o  For  a  relative  length 
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ratio  of  L/b»'9  four  models  were  made:  one  for  each  of  the  following  values 
of  m=b/8,  m-'«0o3j  0»5,  0„7j>  and  0„9»  They  were  constructed  with  1/2  inch 
marine  plywood  and  faced  with  22  gauge  galvanized  sheet  metalo  The  three 
dimensional  models  were  built  with  m  values  of  0„3,  0„p,  0„7,  and  0,9.. 
In  each  "m"  group  two  models  were  constructed  with  L/b-0o25  and  one  model 
with  L/b«0e5o  The  main  construction  wa3  1/2"  marine  plywood-  The  barrel 
was  formed  with  galvanized  sheet  metal,  and  one  side  of  ogse  of  the   L/b«=Oo25 
models  was  faced  with  IuciteD  Figure  6a  illustrates  the  three-dimensional 
bridges,,  Shewn  are  the  four  models  with  L/b=0.25  and  m=0»35  0o5,  0*7,  and 
0o9c  The  book  and  hammer  are  included  to  show  perspective^,  with  this 
combination  of  models  we  were  able  to  test  each  of  the  openings  ra«0o3.c  0=5? 
0o7.«  and  0o9  for  relative  lengths  L/b  of  0,  0=25,  0»50,  0„75>  and  loOOo 

3=  Boundary  Roughness  Analysis 

The  actual  tests  to  determine  the  maximum  backwater ff  were  run 
under  two  different  boundary  roughnesses,,  The  first  roughness^  which  will 
be  called  the  smooth  boundary,  consisted  of  the  steel  walls  of  the  flumee 
/•he  walls  were  finished  with  an  epojrcy  resin  paint e  It  was  determined  that 
the  Siiitioth  flume  produced,  a  Mannings  n  value  of  o01100  This  value  was 
not  representative  of  any  natural  physical  condition^  It  was  d<jcidel  to 
run  a  second  series  of  tests  in  a  boundary  roughness  which  would  simulate 
a  more  natural  condition,,  Assuming  a  scale  of  1/10  between  model  and 
prototypes,  a  Mannings  n  between  0o02  and  0o03  in  the  flume  would  have  been 
desirableo  This  would  correspond  to  field  values  of  approximately  0  03 
to  0o05  respectively,,  The  flume  was  lined  with  a  series  of  1/4  inch 
aluminum  rods«  Two  layers  of  rods  were  placed  on  the  bed  of  the  flume: 
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Figure  5a)  Instrument  Rack 


Figure  5b)  Manometer  Bank 


Figure  6a)  Models 


Figure  6b)  Flume  wifeh  Roughness  Bare 
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a  bottom  layer  of  'ucwgiivtxttml  bars  placed  12  inches  on  canter  and  a  top 
layer  of  transverse  bars  6  in0  on  center0  Along  the  side  wall,  the:.-   ras 
one  layer  of  vertical  bars  6  in<>  on  canter  placed  1/4  in°  from  the 
The  bottom  bars  were  tied  together.  The  vertical  bars  were,   tied  at  ihe  bot- 
tom end  to  the  transverse  bars  and  were  clumped  to  the  wails  above  the  free 
surface,.  This  roughness  pattern  is  shown  in  figure  6bD  After  a  series  of 
normal  depth  test  runs  were  made,,  it  was  found  that  this  particular     srn 
gave  a  Manning's  n  of  0o0236o  It  was  also  found  that  steeper  slopes  and  greater 
depths  could  be  used  without  going  out  of  the  test  rangec 

In  order  to  test  the  consistency  of  the  uniform  flow  depth,,  for  a 
given  slope  and  flow,  comparisons  to  several  well  known  experimental  results 
were  made0  The  Darcy-Weisbaeh  friction  factor  and  the  Reynolds  Number  for 
each  uniform  flow  condition  wa-j  eaculatsdo  For  the  "smooth"  tests*  .:.;  exper- 
imental frict.tcn  factors  were  compared  to  the  theoretical  values  obtained  by 
adapting  the  Biasius  and  Prandtl-Von  Kenna*}  formulas  for  flow  in  smooth  pipes 
to  the  rectangular  open  channel,,  Figure  7  shows  a  plot  of  the  Darcy  ,,'eisbach 
friction  factor  vs0  the  Reynolds  Number  for  both  t!he  smooth  and  rough  test 
data,,  The  rough  data  has  been  broken  down  according  to  constant  flow  lines  and 
constant  hydraulic  radius  lines„ 

Sayre  and  A3±>ertsori*-nave  presented  a  compreshensive  report  of  the 
effect  of  roughness  elements  in  rigid  open  channalsc  They  state  that  a  rough- 
ness parameter  ^  (chi)  which  depends  "on  the  size,  shape  and  spacing  of  the  rough- 
ness  elements" 9   should  completely  describe  the  boundary  roughness*.  The  true 
value  of  yi  depends  on  'whether  or  not  1)  the  boundary  is  hydrodynamically  rough- 
negligible  viscous  effects,  and  2)  the  channel  is  sufficiently  wide  such  that 
any  appreciable  side  wall  effects  are  essentially  eliminated^  The  general 


2U* 


resistance  formula  for  rough  flow  given  by  Sayre  and  Albertson  is 


- 


j±r  =    -OS  \o3J6/x  6) 

According  to  the  method  they  have  described  for  determining  the  value  of  * » 
the  pattern  of  1/4  inch  aluminum  reds  used  in  this  research  gave  a  %   value 
of  o 0126c  The  value  of  6o06  agrees  very  well  with  the  present  worko  Figure 
8  shows  a  graph  of  the  roughness  function(  C/fQ   )  vs0  the  relative  roug] 
ness  y<>/%   for  sorii3  of  the  rough  normal  depth  data0 

Several  velocity  profiles  were  taken  at  a  condition  of  m* 
slope  and  maximum  flow0  iith  the  value  of  jt    and  the  shear  velocity  defined 


as   \/£>/f   -   /i  3°  5    a  diniensionless  velocity  profile  was  drai 
The  equation  deseribiag  this  profile  is 

—£==,   ==  S.OS   ioq  yZ,    +-   4-6  -7} 

J  l"/p  J'°  /7- 

Figure  9  shews  the  graph  of  this  equation  aid  compares  it  to  the  similar 
one  defined  by  Sayre-"*-0  The  difference  in  the  constant  may  be  due  to  the 
intense  wall  effects  which  were  present 0 

;ure  3.0  shows  a  general  resistance  diagram  for  open  charK 
flow*  It  is  similar  in  nature  to  the  famous  Moody  diagram  for  pipe  flowc 
The  curves  that  are  plotted  are  those  suggested  by  Sayreo-^The  smooth  and 
rough  test  values  have  been  plotted  for  comparisonc 

PRELIMINARY  T£STINSi 

Several  preliminary  tests  were  run  in  the  small  6"  by  12»  flumec 
The  results  of  the  two  dimensional  weir  tests  were  put  in  graphical  form,  by 
plotting  the  coefficient  of  discharge  vso  the  contraction  ratio  m"  ith  the 
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GENERAL    RESISTANCE     DIAGRAM    FOR    UNIFORM 


FLOW    IN    OPEN    CHANNELS 

(SAYRE) 
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Froude  Number    Fb  as  a  parameter,,     In  addition,  the  relationship  be? 

Wo   and  jj,/u  was  plotted  in  a  sijrtilar  manner  <>      (i<,e«  m8  as  the  variable   3nd 

Po  as  the  parameter).. 

The  tv/o-djjaensional  case  was  extended  to  the  three-diraonsi ■:. ;  -.1 

case  by  us     emi-circular  arch  bridge  models  of  the  same  b/B  ratio  and 

model  length  L  of  24  inches*  A  compai*ison  of  the  two  and  three-dimensional 

tests  are  compared  in  figure  110  It  is  interesting  to  note  that  at  Froude 

Numbers  less  than  0»5  the  effect  of  length  was  almost  negligible ,  However 
at  higher  Froude  Numbers  the  three-dimensional  tests  exhibited  a  smaller 

value  of  Cd  and  a  larger  backwater  ratio     0 

The  two-dima  semicircular-rough  tests  were  analysed  and  plotted 0 
When  the  corresponding  smooth  tests  were  compared  3  the  differences  were 
found  to  be  negligible e  This  indicated  that  the  boundary  roughness  was  not 
an  influencing  parameter  at  Froude  Numbers  less  than  0„5o  It  was  possible 
that  the  small  scale  effects  due  to  increased  surface  tension  could  result 
in  such  a  misleading  conclusions  It  was  therefore  necessary  to  verify  this 
conclusion  on  a  larger  scale 0 

EXP  5RMfiNTAI,_  PRQgggUgg  i 

When  the  smooth  flume  tests  were  started,  the  procedure  was  to 
arbitrarily  select  a  slope  and  a  flows  and  then  adjust  the  tailgate  until 
a  satisfactory  normal  depth  was  obtained »  This  proved  to  be  a  very  tedious 
and  unsatisfactory  methodo  Before  the  rough  testing  was  begun,  an  effort 
was  made  to  determine  an  exact  relationship  between  the  several  different 
variables  required  to  produce  a  normal  depth  profilec  A  series  of  24  different 
normal  depth  tests  were  run» 
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FIGURE      12 
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From  these  preliminary  tests,  an  average  value  of  Mannings  n  of 
0,0233  was  computed^  A  calibration  chart  for  selecting  normal  depths  was 
made.  This  calibration  chart  is  shown  in  figure  I2C  By  using  these  curves 
to  predetermine  the  slope,  tail-gate-setting  and  normal  depth  to  gi\ 
desired  Froude  Number  at  a  given  discharge,  a  vast  amount  of  time-co tsuraing 
work  was  eliminated.  This  method  was  used  throughout  the  series  of  rough 
tests  and  proved  to  be  very  satisfactory ,  Only  a  few  minor  adjusting  lbs  were 
needed 0 

The  test  procedure  consisted  of  presetting  the  flov.,  slop 
tailgate  and  then  testing  the  various  b/B  and  L/o  models  at  these  normal 
depth  conditionSo  All  depths  were  measured  relative  to  the  flume  bottom,, 
Depths  were  read  along  the  centerline  until,  the  maximum  upstream  point  had 
been  reached  and  passed*  Measurements  were  then  taken  along  the  downstream 
centerline  until  the  minimum  point  had.  been  reached  and  passed*  This  pro- 
cedure was  used  for  95  smooth  tests  and  IK)  rough  tests,  A  more  extensive 
study  of  the  surface  topography  and  j-slocity  profiles  were  performed  on  a 
few  runs, 

T£3TS_AJ©J^ULTS: 

.e   expe  imental  values  of  yjij^for  semi-circular  constrictions  in 
a  smooth  rectangular  channel  were  plotted  vs,  the  contraction  ratio  m8  and 
is  shown  in  figure  13aQ  In  a  similar  manner ,  the  discharge  coeffici 
Cd  for  the  smooth  flume  tests  is  shown  in  figure  13b,  The  equation  used 
to  calculate  this  Cd  is  shown  in  the  figure. 

Since  the  smooth  tests  included  only  the  results  of  the  two  di~ 
raensional  models,  it  was  needed  to  investigate  further  the  effects  of  length 
as  well  as  roughness  in  the  rough  tests.  In  order  to  test  at  low  Froude 
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Numbers  &  very  mild  slope  was  requiredo  Due  to  the  smooth  walls,  ii 

difficult  to  obtain  stable  flow,.  The  expanding  flow  at  the  downstream  aide 

of  the  constriction  was  often  unstable 0  It  would  deflect  to  one  sic  a  or  the 

other  and  would  seldom  remain  evenly  distributed,,  These  testing  dii  ::':•. culties 

were  probably  the  cause  for  the  scatter  of  smooth  flow  data  points  on  the 

curve  of  the  friction  factor  vs»  the  Reynolds  Number  in  figure  7o 

^°  Rou^h  Boundary  Tests 

The  table  below  shows  the  conditions  which  were  tested  in  the  rough 

channel.  The  X'siadicate  the  desired  normal  depth  conditions  in  which  the 

following  values  of  m  and  L/b  ratios  were  tested: 

ir.  -  b/B  -  0.3,  0c5s  0c?3  0»9 
L/b  -  0,  0o5s  IcO 


?he  experimental  conditions  were  obtained  from  the  calibration  chart  of 
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The  particular  measurements  that  were  taken  on  each  of  the  ?.bove 
mentioned  tests  were  those  required  to  calculate  the  following  quantities: 
The  hydraulic  radius,  the  Reynolds  tfo„  f£0  ,  the  Froude  No*  fF0  ,   the  fric- 
tion coefficient  f*  the  contraction  ratio  m*,  the  discharge  coefficient  Cu, 
the  backwater  ratio  urn  s   the  backwater  superelevation  hi,  the  surface  pro= 
file  ratio  hf/Ah  *  the  length  to  the  maximum  backwater  elevation  L±,   the 
length  to  the  point  of  minimum  depth  L3,  Li  +  L3,  and  Manning's  n<> 
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In  view  of  the  large  amount  of  data  that  was  to  be  analyzed  and 
the  repetitive  character  of  the  calculations,  a  program  was  prepared  for 
processing  the  data  on  the  Royal  McBee  LGP-30  digital  computer 0  The  total 
computer  time  required  to  calculate  all  of  the  above  mentioned  quantities 
for  the  ISO  rough  test  runs  was  a  little  under  nine  hours 0  'without 
computer,  the  time  required  for  processing  the  data  would  have  been  prohi- 
bitive c 

The  backwater  ratio  u  /y  is  shown  in  figure  14  for  the  sei 
circular  constrictions  (L/bO)  in  the  rough  channel*  This  plot  is  similar 
to  that  shown  for  the  smooth  channel  in  figure  13a „  The  graph  cleai 
shows  that  as  mu  approached  unity  y  /y  goes  to  one0  Also  as  m'  ioes  to 
aero  the  backwater  ratio  approaches  infinity <>  The  actual  test  values  are 
not  shown  since  the  curves  of  constant  Froude  Numbers  have  been  graphically 
interpolatedo  The  amount  of  error  produced  during  the  interpolating  process 
was  found  in  most  cases  to  be  less  than  one  per  cento  Similar  plots  have  been 
made  for  the  relative  lengths  L/b  of  0*5  and  lo0o 

A  comparison  of  all  of  the  L/b  plots  for  backwater  ratio  '*,ave   been 
made  in  figure  l6a0  These  curves  were  produced  by  taking  cross-section  at 
m"  values  of  0*3,  0<,5,  and  0„7  from  the  plots  of  y  j\jo   vs.  me0  It  appears 
that  the  Froude  JJo0  and  the  contraction  ratio  are  the  governing  parameters,, 
.Especially  at  lower  Froude  Numbers  (below  0*5) ,  the  influence  of  the  bridge 
length  seemed  to  be  smallo  The  effect  seemed  to  increase  with  a  decrease 
in  the  contraction  ratioc  In  the  case  of  a  small  m89  the  physical  propor- 
tions of  the  constriction  are  closer  to  those  of  a  culvert  rather  than  a 
bridge  opening 0 

Figure  15  shows  the  graph  of  the  discharge  coefficient  vs„  the 
contraction  ratio  with  the  Froude  Number  Fo   as  a  parameter  for  the 
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semi-circular  rough  tests.  These  curves  were  also  interpolated  for 
stant  Froude  Number  linesc  Similar  graphs  were  drawn  for  the  tests 
formed  on  the  other  three-dimensional  models.  The  huT,p  that  appears 
the  Froude  ilumher  3.ines  of  0,25  to  0,60  was  a  phenomena  which  appeared  in 
all  of  the  plots  of  the  rough  tests*  A  comparison  of  the  several  le 
ratios  was  also  done  by  taking  cross- section  at  constant  m9  values,  A  ty- 
pical cross-section  at  m!  ■  0,7  is  shown  in  figure  16b,  This  graph  as  well 
as  other  similar  ones  strongly  reveals  the  fact  that  the  bridge  length  is 
relatively  unimportant  and  can  for  all  practical  purposes  be  disregarded* 

In  figure  17  the  results  of  both  smooth  and  rough  tests  are  com- 
pared by  the  method  of  cross-sections „  These  curves  verify  the  conclusion 
.made  from  the  small  flume  tests,  that  below  a  Froude  Number  of  0o5  ': 
backwater  produced  by   a  given  constriction  is  essentially  the  same 
smooth  and  rough  boundaries0 

In  order  to  completely  describe  the  centerline  profile  it  is 
desirable  to  have  an  estimate  of  the  distance  from  the  upstream  face  of  the 
constriction  to  the  point  of  maximum  backwater  elevation.  This  distance 
is  referred  to  as  Lt_0  Because  of  the  flatness  of  the  surface  profil 
the  vicinity  of  the  maximum  point,  it  was  extremely  difficult  to  get  an 
exact  measurement  of  Li,  The  actual  measurements  taken  could  have  b 
in  error  by  as  much  as  +  0,5  feet.  However 3   with  the  large  amount  of  data 
which  was  available,  it  was  possible  to  study  Lj,  on  an  average  basis.   Aver- 
age values  of  L\  were  calculated  for  several  combinations  of  b/B,  L/1  ,  L/B, 
etCo  In  this  manner,,  it  appeared  that  the  variable  bridge  length  and  the 
change  in  m*  were  of  the  same  order  of  magnitude  as  the  experimental  error. 
The  most  consistent  relationship  was  foiand  by  plotting  the  dimensioniess 
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ratio  I/j/b  vsc  the  Froude  No»  with  m=b/B  as  the  parameter „  This  relation- 
ship is  shown  in  figure  13a„  The  valuefflcf  Li  obtained  from  the  smooth  tests 

also  compared  favorably  with  figure  l?ia0  In  a  similar  manner  it  was  found 
that  the  length  L3/L.3  (distance  from  the  maximum  point  to  the  minimum  point) 
varied  only  with  the  constriction  geometryc  The  average  values  of  Lp-I/j/b 
are  plotted  vs0  m  -  b/B  with  L/b  as  a  parameter  in  figure  18bc 

Several  other  investigators  have  used  the  ifroude  Number  at  sec- 
tion 3  (see  figure.!)  as  an  estimator  of  the  maximum  backwater a     Others  have 
used  W$     as  a  controllj-ng  parameter  in  making  indirect  measurement"  of 

od  discharges,.  Due  to  the  extremely  irregular  flow  pattern  at  the  mini- 
mum point  it  would  seem  that  the  use  of   iF3  may  be  very  misleading „  In 
the  present  research,  the  normal  depth  Froude  No0  Wo     was  found  to  be  s. 
very  reliable  estimator  of  y  / u    0  In  order  to  test  the  variability  of  TFg 
with  fFb  a  correlation  curve  of  fl^/ii^vso  W&     was  preparedc.  Tliis  curve 
is  shown  in  figure  190  Below  a  Froude  Number  of  0o5  the  correlation  was 
goodo  However 3   above  ffe,  a0o5  the  depth  V3  wa3  often  below  the  critical 
depth  and  the  correlation  of  ff^/fife t0  ^  wa3  very  poor0  The  scatter  seemed 

to  incx'ease  with  increasing  values  of  L/bu  Therefore  only  the  test  results 
of  the  L/b«0  tests  are  shown*  If  used  with  caution  these  curves  can  be  used 

to  estimate  the  minimum  depth  y*0     It  appears  from  this  curve  that  his  a 
much  more  reliable  a^ti^ator  than  fi^  o 
Go  Segment  Analysis  ana  CcBtperlsions 

With  the  introduction  of  mss  the  assumption  was  made  that  if  pro- 
perly interpreted  the  backwater  produced  by  constrictions  of  the  same  m" 
would  be  equal  regardless  of  the  physical  geometry  of  the  actual  constric- 
tiono  In  order  to  varify  this  assumption  test  data  on  constriction  geome- 
tries other  than  a  semi-circle  was  neededo  A  series  of  90  tests  were  run 
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in  the  small  preliminary  flume  on  two-dimensional  segment  weirs  witl 
fl~d/r  value  of  0o5o  (see  figure?  )„  The  data  obtained  were  reanal, 
in.  terms  of  m«  0  These  tests  were  run  in  the  rough  channel  which  hac 
Mannings  n  of  0o0201o  After  the  results  had  been  plotted  in  the  fox-a  of 
y  /y    vso  m"  with  W0    as  a  parameter,  a  ©oaEpscriLeea  was  raade  with  tho  two- 
dimensional  rough  tests  in  the  large  flume.  The  results  of  the  comparison 
we.2'€;  very  goodo  Inspite  of  the  fact  that  each  set  of  cur  nter- 

polated_,  the  small  differences  in  the  plots  could  ee.s: .'     .ttributod  to 
experimental  and  grapljical  errors,. 

In  a  similar  manner,  the  vertical  board  data  given  by  Liu^'  was 
reanalyzed  to  fit  the  plot  of  l>  ju   vs»  m*0  These  tests  were  run  in  a  wider 
flume  with  a  different  roughness  pattern0  Their  roughness  produced  a  Man- 
ning's n  of  Co024o  The  results  were   compared  to  the  semi-circular  data 
and  the  segment  data0  Again  the  differences  were  extremely  small  an:!  con- 
tributable  to  experimental  error0  The  authors  feel  that  it  is  extremely 
interesting  to  note  that  the  test  data  taken  by  three  investigators  in 
three  different  flumes  and  under  three  completely  different  set  ups  pr 
duced  almost  identical  resultSo  This  clearly  verifies  that  as  defined  the 
contraction  ratio  m!  is  essentially  an  ail  inclusive  term0  Of  course  the 
data  compared  were  those  where  the  eccentricity  was  zero,  the  skew  was  aero, 
and  the  entrance  was  sharp „  It  is  still  necessary  to  apply  correction  terns 
for  these  conditions* 

It  would  seem  that  due  to  the  similar  results  mentioned  above 
there  should  be  some  relationship  between  the  backwater  ratio  t/,/yd  ,  the 
Froude  Noc  ?F0     ,   and  the  contraction  ratio  ml  This  relationship  should 
be  applicable  to  all  constriction  geometries 0  As  mentioned  previously  in 
the  analysis,  a  similarity  was  noticed  between  the  several  different 
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backwater  equations,,  The  term  (fFb/^'j  appeared  in  all  of  the  solut: 
of  'j  /%  °     &1  general  it  appeared  that 


M 


y0 


-  c 


(  m'J 


2/3 


(30) 


where  C  is  a  coefficient  which  would  take  in  the  effects  of  the  disc:  arge 
coefficient,  approach  velocities,  nonuniform  velocity  distributions  and 
other  empirically  determined  factors,-  Equation(30)  is  actually  the  equa- 
tion of  a  straight  line  on  logarithmic  paper  with  a  slope  of  5o  A  t 
of  $0  semi-circular  L/b*0  test  values,  44  vertical  board  values  ( Cole rado  ) 
and  50  segment  values  were  plotted  in  the  form  of  y,/y,-' vs0  (ifo/m)  anc  are 
shown  in  figure  20o  The  value  of  tyy  - 1  was  used  in  order  to  expand  the 
scale  of  the  backwater  ratic  It  is  quite  apparent  that  the  data  collapsed 
into  one  general  straight  line  relationship a 

The  method  of  least  squares  was  applied  to  a  random  sample  of  the 
144  test  points  to  determine  the  straight  line  relationship0  After  solving 
for  t      and  Cv  eq0  (30)  became 
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jtSqo  (31)  is  a  very  simple  and  easy  solution  for  the  backwa 
produced  ay  any  type  of  constriction,,  In  actual  practice,  this  equation 
will  give  as  good  an  estimate  of  the  maximum  backwater  y^  as  any  previous- 
ly suggested  method. 
Do  Surface  Topography  and  Velocity  Diagrams 

In  order  to  complete  the  analysis  of  the  maximum  backwater,  addi- 
tional studies  were  made  of  the  velocity  distributions  and  the  surface  pro- 
files., These  studies  were  made  for  the  condition  of  a  sharp  crested  semi- 
circular  constriction  with  m=0<3  at  a  Froude  Noc  ff0   of  0o5° 
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A  detail  of  the  surface  topography  both  upstream  and  downstream 
was  observed.   The  result  of  this  study  is  shown  in  figure  21,.  The  u  bens 
shown  are  the  depths  in  centimeters „  Only  a  detail  of  one  half  of  the  sv 
face  is  shown  since  the  other  side  is  essentially  symetric0  The  graph  shows 
lines  of  equal  surface  elevation,,  The      "'ine  surface  profile  ic  shown 
on  the  top  =  It  is  interesting  to  note  that  the  actual  maximum  water  surface 
elevation  is  not  along  the  centerline.,  but  on  the  upstream  face  of  the  abut- 
ment o  Al      this  may  be  ejected,  the  actual  magnitude  of  the  differ- 
ence in  elevation  between  the  d     maximum  elevation  and  the  actual  maximum 
is  the  important  question o  The  actual  maximum  shoreline  elevation  was  found 
to  exceed  the  maximum  centerline  elevation  by  as  much  as  5$  of  the  center line 
depth.  This  fact  was  verified  in  the  surface  topographies  taken  at 
conditions,,  Liu",  as  well  as  Herbich^  gave  similar  surface  topographies  of 
other  geometry  constrictions,  and  the  difference  in  water  surface  elev  tions 
was  again  found  to  be  about  5k   of  the  centerline  depth,,  In  general,  ic  seems 
that  any  e       of  the  maximum  centerline  depth  yn_  should  be  increased  by 
y,i   to  get  the  maximum  shoreline  elevation.. 

In  addition  to  the  surface  topography,  several  velocity  prciiles 
were  taken*  Traverses  were  taken  with  the  Prandtl  Tube  at  four  sections 
with  the  model  mO„3  and  L/b=0  and  a  Froude  Mo.  Fib  of  0o5o  The  first  section 
was  in  to  backwater  region  of  essentially  uniform  flow0  The  second  i  as  at 
the  section  of  maximum  backwater,  the  third  at  the  vena  contracta  an<  the 
fourth  at  the  section  of  minimum  depth.  At  each  section  a  vertical  •=  elocity 
traverse  was  taksn  at  1  £t„,  2  £t„  and  2C35  ft«  both  left  and  right  of  the 
centerlineo  At  the  vena  contracta  they  were  taken  at  the  <h    ,  0„5  feet  and 
1  foot  left  and  righto  In  general^  a  more  detailed  traverse  was  taken  at 


the  centerlinec  From  these  measurements ,  plots  of  equal  velocity  line3  ware 
prepared  for  the  several  sections,,  A  composite  picture  of  the  isovel  dia- 
grams is  shown  in  figure  22c  Only  one  half  of  the  diagrams  are  shown  due  to 
symetry.  All  of  the  diagrams  were  integrated  by  a  planimeter  and  the  discharge 
was  checked  against  the  ven turd-meter  discharge.,  They  all  checked  within  l#o 
The  kinetic  energy  coefficient   ex  and  the  momentum  coefficients  p     were 
calculated  for  the  section  of  uniform  flow  and  were  found  to  be  oi'=\  '.-■    and 
8'  -    1.18  respectively o 

Aa  estimate  was  also  made  of  the  force  required  by  the  bridge  to 
produce  the  resulting  backward 0  This  way  done  by  applying  the  momentum  equa~ 
tion  in  the  integral  form  between  the  section  of  maximum  backwater  and  the 
vena  contraeta0  By  integrating  the  isovel  diagrams  of  figure  21  and  apply- 
ing the  momentum  eouation,  the  required  bridge  force  was  found  to  be  4c  54  lbs, 
if  a  similar  calculation  is  made  on  a  similar  prototype  bridge  with  a  model- 
prototype  scale  t'i*  1/2CV  the:  tr-ldge  fcr«s  wcalci  J»b  ??6,000  lbs0 

CONCLUSIONS  AND  jjEgIGN  j^CQI^NDATIQNS : 

The  most  important  variables  in  determining  the  maximum  backwater 
are  the  normal  depth  Frcude  Number  fF0    and  the  contraction  rat-io  mK     As 
defined  the  contraction  ratio  can  be  used  for  any  and  all  types  of  bridge 
con strict! on So  The  boundary  roughness  as  well  as  the  bridge  length  for  Froude 
Number   RF0   less  than  0o5  are  relatively  unimportant,  and  their  effects 
for  all  practical  purposes  can  be  neglected.,  The  best  approximation  co  the 
backwater  ratio  for  semi-circular  arch  bridges  is  given  in  figure  14c  Equa= 
tion(21)  can  be  used  to  calculate  y Jy    by  obtaining  the  discharge  coeffi- 
cients from  figure  15-,  A  more  practical  first  approximation  k.  the  maximum 
backwater  is  given  by  the  curve  of  figure  20  or  equation  (31) „  It  is 
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5io 
recommended  that  5$  of  this  maximum  depth  be  added  to  the  corresponc 
center line  elevation  to  get  the  maximum  shoreline  elevationc  This  eleva- 
tion should  occur  in  the  vicinity  of  the  bridge  embankment o 

In  determining  the  backwater  produced  by  a  new  single  spar,  semi- 
circular symetric  arch  bridge  where  the  spriagline  of  the  arch  is  at  the 
bed  of  the  stream,  the  following  design  procedure  is  recommended* 

lo)  Plot  the  normal  depth  on.  a  section  view  of  the 

stream  eross-section  where  the  bridge  is  to  be 
built o 

20)  Superimpose  the  proposed  bridge  design  on  this 

section  view0 

3o)  Determine  the  value  of  m  »  b/Bo 

Calculate  y  /r  and  get  the  value  of  Cm  from 
figure  3  for  the  curve  6jr=0o     When  the  center 
of  curvature  is  below  the  spriagline^,  calcu- 
late d/r  and  us®  the  respective  curve  to  obtain 
Cm0 

5<0  Calculate  the  norma],  depth  Froude  Number-  FFS  „ 
(The  discharge  should  be  given  and  the  average 
velocity  V0    can.  be  determined  from  the  con- 
tinuity equation.,) 

6C)  Gaiculate  m"   »  ££,mo     (The  value  of  si*   could  be 

checked  by  planiiaetering  the  areas  and  getting 
the  ratio  AjA^  directly  <> ) 

7o)  With  m'  and  W0     get  the  value  of  y ly   from  fi= 
gure  14 o  A  more  approximate  value  can  be  ob- 
tained from  figure  20  or  equation  (31) = 


Calculate  yx  arid. 

shore 

..ith     FF0  L]_« 

3  15b  .-. 

a 

- 

If  an 

■La.   the  oc 

ilated0     The  following  procedure  . 

3d; 

i 

.scharge,  calculate  the  no:         of 
flow  - 

2,,'  ting  surface  ele- 

imura  centerline  depth  y]  can  be     -  sd 

iing  the  difference  between     saax„  elevation 

eiev.  of  the  bed  of 

Note:  This  accounts  for  ths  y£   difference  in 
depth  between  tl     terline  and  the  stre 
banks0  This  step  ;an  be  disregaxded  as  a  safety 
f  ac      estima     '  he  minimum  bridge  ar 

3o;  Calculate  the  normal  depth  Froude  Number  (t0 

4o/        lu     and.  ffc    obtain  the  value  of  the  con- 

traction  ratio  m°  from.  figure  14  or  20  &;.'  by 

oqu^tlon  (31°) 
5>/  Calculate  the  total  normal  depth  flow  area  A- 


6:)  w.'.tn  m1  and  Aq  d   calculate  the  desir  .         imum 

.e  value 
ratio  ia'     ;  corre 
determine 
Figure  3 

sc      to  the  radius 

~:. 

If  i  ..  re  used  as  a  method  for 

.teg  of  flood  discha- 
recc;  tended  ,  It  is  is  a 

a  direct 

i  the  3».?diLum 
bae. 
20)  depth  y-,  c 

3o> 

the  bi  aection  vi 

;alcula;.  normal  c  x&ther 

equation  11  areas 

gett    ■  bly) 

its  the  Froude  Number  from  equation  31° 
6o)     With     ffv>     and     ya    use  equation.  19  to  get  the 

first  evStiniate  of  the  di3charge0 
7o)     With     Wo    and  m°   get  the  value  of  the  discharge 
coefficient  from  figure  15 o      (Equation  20  can 
also  be  usedD) 
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^  ical  depth  flow  area  of 

ictioric. 

m  or  fli 

si: 

b  Width  of  I 

C  /T  roughness 

Cd  ioient  of  rge 

Cm  stion  coefficient o 

d  fr«« 

f 

g 

Ah 


1  L  jfch  of  the  bridge  paralle_  to 

the  direction  of  flow* 

**1  L  The  perpendicular  distance  f 

the  upstream  face  of  the  bri 
to  the  maximum  backwater  elevationc 

^3  L  The  perpendicular  distance  from 

the  upstream  face  of  the  bridge 
to  the  minimum  surface  elevationo 


bure  of 

Numb  3'. 

r  at 

3  minim 

i/eisba< 

tion 

Accele 

ration 

of 

sncs  betwe< 

sn  i 

:.he  max? 
elevat; 

580 

M  An  infinite  series  of  powers  of 

the  maximum  depth  to  radius  ratioo 

M^  Open  channel^  mild  slope      iter 


m  zh  ratio  b/Bo 

ra1  Channel  contraction  ratio0 

n  L'  Mannings  roughness  coef f ici  - 

0  L3-T1  Total  flow, 

q  1?   "■*-  Portion  of  the  total  flov 

could  pass  through  the  bri 
without  contractiono 

R  L  Hydraulic  radius c 

g>  Reynolds  Number  *  vy 

RQ  Reynolds  dumber  =»  ^^/i? 

r  L  Radius  of  the  arch. 

y0  L  Depth  of  normal  unconstri     Iow0 

y\  L  Depth  of  the  maximum  back 

j 2  Z>  Depth  at  the  vena  contracta: 

7^  L  Depth  at  the  minimum  surf ac , 

elevationo 

Symbol  for  the  ratio  ya/r  « 

ex'  Kinetic  energy  coefficient: 

B  Symbol  for  the  ratio  d/rD 

B'  Momentum  coefficient o 

-f>  L"T  Kinematic  viscosity  of  the  fluido 

o  FT  V*  Fluid  mass  density,, 

V0  FL  Shear  intensity  acting  on  the 

channel  bedc 

?C  L  Roughness  parameter  suggested  by 

Sayre0 


